#include <math.h>
#include <stdio.h>

//#define PI 3.14159265359

float sx, sy;

float sdCircle(float px, float py, float r) {
	float dx = px - sx, dy = py - sy;
	return sqrtf(dx * dx + dy * dy) - r;
}

float opUnion(float d1, float d2) {
	return d1 < d2 ? d1 : d2;
}

#define T px + scale * r * cosf(theta), py + scale * r * sin(theta)

float f(float px, float py, float theta, float scale, int n) {
	float d = 0.0f;
	for (float r = 0.0f; r < 0.8f; r += 0.02f)
		d = opUnion(d, sdCircle(T, 0.05f * scale * (0.95f - r)));

	if (n > 0)
		for (int t = -1; t <= 1; t += 2) {
			float tt = theta + t * 1.8f;
			float ss = scale * 0.9f;
			for (float r = 0.2f; r < 0.8f; r += 0.1f) {
				d = opUnion(d, f(T, tt, ss * 0.5f, n - 1));
				ss *= 0.8f;
			}
		}

	return d;
}

int ribbon() {
	float x = (fmodf(sy, 0.1f) / 0.1f - 0.5f) * 0.5f;
	return sx >= x - 0.05f && sx <= x + 0.05f;
}

//void tree1() {
//	int n = 3;
//	for (int j = 1; j <= n; j++) {
//		int s = 1 << j, k = (1 << n) - s, x;
//		for (int y = s - j; y >= 0; y--, putchar('\n')) {
//			for (x = 0; x < y + k; x++)
//				printf("  ");
//			for (x = 0; x + y < s; x++)
//				printf("%c ", '!' ^ y & x);
//			for (x = 1; x + y < s; x++)
//				printf("%c ", '!' ^ y & (s - y - x - 1));
//		}
//	}
//}

void tree2() {
	int n = 2;
	for (sy = 0.8f; sy > 0.0f; sy -= 0.02f, putchar('\n'))
		for (sx = -0.35f; sx < 0.35f; sx += 0.01f)
			putchar(f(0, 0, PI * 0.5f, 1.0f, n) < 0 ? '*' : ' ');
}

void tree3() {
	int n = 2;
	float zoom = 1;//argc > 2 ? atof(argv[2]) : 1.0f;
	for (sy = 0.8f; sy > 0.0f; sy -= 0.02f / zoom, putchar('\n'))
		for (sx = -0.35f; sx < 0.35f; sx += 0.01f / zoom) {
			if (f(0, 0, PI * 0.5f, 1.0f, n) < 0.0f) {
				if (sy < 0.1f)
					putchar('.');
				else {
					if (ribbon())
						putchar('=');
					else
						putchar(
								"............................#j&o"[rand() % 32]);
				}
			} else
				putchar(' ');
		}
}
